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00 ■ Abstract. The study of polynomial solutions to the classical Lame equation 

1 in its algebraic form, or equivalently, of double-periodic solutions of its Weier- 

' strass form has a long history. Such solutions appear at integer values of the 

s ! , spectral parameter and their respective eigenvalues serve as the ends of bands 

, in the boundary value problem for the corresponding Schrodinger equation 

■ with finite gap potential given by the Weierstrass p-function on the real line. 

In this paper we establish several natural (and equivalent) formulas in terms of 
hypergeometric and elliptic type integrals for the density of the appropriately 
scaled asymptotic distribution of these eigenvalues when the integer-valued 
spectral parameter tends to infinity. We also show that this density satisfies a 
Heun differential equation with four singularities. 

ji ■ 

■M ■ 1. Introduction and Main Results 

C3 • 

The algebraic form of the classical Lame equation [STJ Chap. 23] was introduced 
i—i i by Lame in the 1830's in connection with the separation of variables in the Laplace 

| equation by means of elliptic coordinates in K'. Lame's equation is given by 

' . .... 

where Qi(z) is a real degree I polynomial with all real and distinct roots and V(z) 

is a polynomial of degree at most I — 2 whose choice depends on the type of solution 

to {1} one is looking for. In the second half of the nineteenth century several famous 

mathematicians including Bocher, Heine, Klein and Stieltjes studied the number 

and various properties of the so-called Lame solutions of the first kind (of given 

| degree and type) to equation (fTJ) . These are also known as Lame polynomials of a 

certain type. Such solutions exist for special choices of V(z) and are characterized 

by the property that their logarithmic derivative is a rational function. For a given 

Q(z) of degree I > 2 with simple roots there exist 2 l different possibilities for Lame 

. £^ | solutions depending on whether these solutions are smooth at a given root of Q{z) 

^ ■ or have a square root singularity, see [HI E] for more details. 

A generalized Lame equation [21] is a second order differential equation of the 

form 

d 2 , . d 



C3 



Q{z) dz^ +P{z) d~z +V{z) j S{z)= °> (2) 
where Q{z) is a complex polynomial of degree / and P{z) is a complex polynomial 
of degree at most l — l. As it was first shown by Heine [TU] for a generic equation of 
either of the forms JT]) or ^ and an arbitrary positive integer n there are exactly 
Cn 2 ) polynomials V(z) such that S(z) is a polynomial of degree n. Below we 
concentrate on the most classical case of a cubic polynomial Q(z) with all real 
roots as treated in e.g. [21]. Already Lame and Liouville knew that in this case the 
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(unique) root of V{z) is real and located between the minimal and the maximal 
roots of Q{z). The latter property was further generalized by Heine, Van Vleck and 
other authors. Note that if degQ(z) = 3 then V{z) is at most linear and that for 
a given value of the positive integer n there are at most n + 1 such polynomials. 

Before formulating our main results let us briefly review - following mostly [131 
§3] and [14] - some necessary background on the version of the Lame equation 
used below and its spectral polynomials. Setting, as one traditionally does, Q{z) = 
4(z — e±)(z — e-2){z — e^) with e± > €2 > e$ we can rewrite equation (fTJ) as 



(The chosen representation for the linear polynomial V(z) will shortly become 
clear.) Notice that several equivalent forms of the Lame equation are classically 
known. Among those one should mention two algebraic forms, the Jacobian form 
and the Weierstrassian form, respectively, see [2"Tl §23.4] and Remark Q] below. 
Equation ^ presents the most commonly used real algebraic form of the Lame 
equation which is smooth while the other (real) algebraic form is sing ularQ 

A Lame solution of the first kind to equation ([3]) is a solution of the form S(z) = 
(z — ei) Kl (z — e-i) K ^{z — e 3 ) K3 S(z), where each m is either or | and S(z) is a 
polynomial. Lame solutions which are pure polynomials - i.e., for which k\ = K2 = 
«3 = - are said to be of type 1, those with a single square root are of type 2, those 
with two square roots are of type 3 and, finally, those involving three square roots 
are said to be of type 4- ("Types" are sometimes called "species", see, e.g., [2"T].) 
One can easily check that ([3]) has a Lame solution if and only if n is a nonnegative 
integer. Moreover, if n is even then only solutions of types 1 and 3 exist. Namely, 
for appropriate choices of the energy constant E one gets exactly distinct 
independent solutions of type 1 and independent solutions of type 3. If n is 
odd then only solutions of types 2 and 4 exist. In this case, for appropriate choices 
of the energy constant E one gets exactly 3 ("+ 1 ) distinct independent solutions of 
type 2 and independent solutions of type 4. In both cases the total number 
of Lame solutions equals In + 1, which coincides with the number of independent 
spherical harmonics of order n. 

Let R n {E) = rijUo 1 ^ - Ej), n eZ+, denote the monic polynomial of degree 
2n + l whose roots are exactly the values of the energy E at which equation Q has 
a Lame solution. These polynomials are often referred to as spectral polynomials 
in the literature; their study goes back to Hermite and Halphen. The most recent 
results in this direction can be found in [T31 E3 E3 [2U] and [8, 9J. In particular, 
[13j contains an excellent survey of this topic as well as a comprehensive table with 
these polynomials (and their modified versions) correcting several mistakes that 
occurred in previous publications. Article [9] is apparently the first attempt to 
give a (somewhat) closed formula for R n (E); for this the authors use yet another 
family of polynomials which they call elliptic Bernoulli polynomials. Since explicit 
formulas for R n (E) seem to be rather complicated and difficult to handle, in this 
paper we study the asymptotics of the root distribution of appropriately scaled 
versions of R n (E) (Corollary [1]). In spite of a more than 150 years long history of 
the Lame equation the only source discussing questions similar to ours that we were 
able to locate is [B]. We should also mention that the results below are actually 
much more precise and therefore supersede those of loc. cit. 



The authors thank the anonymous referee for this observation. 
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Let us now introduce eight spectral polynomials i?^ 1,K2:K3 (E) related to the eight 
types of Lame solutions mentioned above, namely 



jy , f R^ Q - {E)Rl h0 (E)Rt°'HE)R n hi (E) when n is even, 

Rn\E) — < l n n nin n n i ill W 

{Rf ' {E)R°r'°{E)R°n ' 2 (E)R r r 2 ' 2 (E) when n is odd. 




Note for example that if n is even then Rfy°>°(E) is the (unique) monic polynomial 
of degree whose roots are precisely the values of E for which equation ^ has 
a pure polynomial solution. 

The results of the present paper actually hold not only for equation ([3]) but also 
for generalized Lame equations (cf. @) of the form 

S(z) = 0, (5) 

where cnj > 0, i = 1, 2, 3, and V(z) is an undetermined linear polynomial. This case 
was thoroughly treated by Stieltjes in [15] . In particular, he proved that for any 
positive integer n there exist exactly n + 1 polynomials V n- j(z), 1 < j < n + 1, such 
that ((5|) has a polynomial solution S(z) of degree n. Moreover, the unique root t n j 
°f Vn.j{z) lies in the interval {e^,e\) and these n + 1 roots are pairwise distinct. 
Consider now the polynomial 

n+l 

i=i 

which is the scaled version of the spectral polynomial R^ 0,0 (E). More exactly, 
for any even n one has that Sp„ (t) — — — " t+z ■ Assume further that the set 

{tn,j : 1 < j < n + 1} is ordered so that t n ^\ < t n ^, < . . . < i n! „+i and associate to 
it the finite measure 



/i„ = ——J2S{z-t n ,,, 



where S(z — a) is the Dirac measure supported at a. The measure /j, n thus obtained 
is clearly a real probability measure that one usually refers to as the root-counting 
measure of Sp n (t). 

Below we shall make use of some well-known notions from the theory of special 
functions that can be found in e.g. pQ. In fact several of our formulas rely on various 
integral representations and transformation properties of the Gauss hypergeometric 
series 

. . , , L(c) T(a + m)T(b + m) z m 

F(a,b,c;z)= 2 F 1 (a,b,c;z)- X 



T(a)T(b) ^ r(c + m) ml' 

where T denotes Euler's Gamma-function, as well as the complete elliptic integral 
of the first kind 

K(C) = / dt =■ (6) 

Recall from [TJ §15.3.1] that if 5R(c) > > then a convenient way of rewriting 
F(a, b, c; z) is given by 

F(a,b,c;z)= ^ (C) f 1 _ t) c-6-i (1 _ u) -*^ (7) 

r(6)r(c - b) j 

The above integral represents a one-valued analytic function in the z-plane cut 
along the real axis from 1 to oo and thus it gives the analytic continuation of 
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F(a,b,c; z). To simplify our formulas let us also define a real- valued function of 
three real variables (xi, x 2 , x 3 ) u>{x\,x 2 ,x 3 ) by setting 

u (x 1 ,x 2 ,x 3 = - — r— j — 8 

\xi\\x 2 - x 3 \ + \x 2 \\xi - x 3 \ 

Note that lu(xi, x 2 , x 3 ) is symmetric in its first two arguments and that 

< oj{x\, x 2 , x 3 ) < 1 whenever x\ < 0, x 2 > and x\ < x 3 < x 2 . (9) 

We should emphasize that the latter conditions on x\, x 2 , x 3 - hence also the bounds 
for cj (xi, X2, x 3 ) given in ^ - always hold in the present context. 
We are now ready to state our main results. 

Theorem 1. For any choice of real numbers e\ > e 2 > e 3 and positive numbers 
a\,a 2 ,a 3 the sequence of measures {/J, n }nez + strongly converges to the probability 
measure /iq supported on the interval [e 3l ei] with density pq given by any of the 
following equivalent expressions: 

(i) 



pqO) 




2 



1 + yl - u{e x - e 2 , e 3 - e 2 , s - 
- e 3 )\s - e 2 \u(ei - e 2 , e 3 - e 2 ,s - e 2 ) 




2y/l - uj{ei - e 2 ,e 3 -e 2 ,s- e 2 ) 2 



1 - y/l - u{ei - e 2 , e 3 - e 2 , s - e 2 ) 2 
where uj is the function defined in ([S]) . 




1 + y_j - ^(ei - e 2 , e 3 -e 2 ,s- e 2 ) 2 
(ei - e 3 )|s - e 2 |cj(ei - e 2 , e 3 - e 2 ,s - e 2 ) 




2-y/l - ujjex - e 2 , e 3 - e 2 , s - e 2 ) 2 
1 - y/l - cj(ei - e 2 , e 3 - e 2 , s - e 2 ) 2 



(iii) 
(iv) 



w(ei - e 2 ,e 3 - e 2 ,s - e 2 ) 



2(ei - e 3 )|s - e 2 |(l + w(ei - e 2 , e 3 - e 2 , s - e 2 )) 



Pl'l 1 1 . 1 - ^(ei - e 2 , e 3 - e 2 , s - e 2 ) 
2' 2' ' 1 + uj(ei - e 2 ,e 3 - e 2 , s - e 2 ) 



^2 [(ei + e 3 - 2e 2 )(s - e 2 ) + 2(ei - e 2 )(e 2 - e 3 ) + (e x - e 3 )\s ~ e 2 \ 



p.'l 1 ^ (gi + e 3 - 2e 2 )(s - e 2 ) + 2(ei - e 2 )(e 2 - e 3 ) - (ei - e 3 )|s - e 2 \ 
" \2'2' ' (e 1 +e 3 -2e 2 )(s-e 2 ) + 2(e 1 -e 2 )(e 2 -e 3 ) + (e 1 -e 3 )|s-e 2 | 

i/iai is, 

1 /" ei da; 

Pq( s ) = TTL / 7= = = = = when e 3 <s< e 2 , 

ZlT Je 2 ^{ei- x)(x - e 2 )(x - e 3 )(x - s) 

1 f 62 dx 

PQ( S ) = / 77 77 w u ^=f ^ en 62 < s < ei - 

^Jea v(ei - a;)(e 2 - x)(a; - e 3 )(s - x) 
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Theorem 2. The density function pq defined in Theorem]^ satisfies the following 
Heun differential equation 

8Q(s)p'^(s) + 8Q'(s)p' Q (s) + Q"(s)p Q (s) = 0, (10) 

where Q(s) = (s — ei)(s — e.2){s — e^,). Both indices of this equation at the finite 
regular singularities e\ , e-i , vanish while its indices at oo equal ^ and | . 

Corollary 1. The root-counting measures for each of the eight normalized spectral 
polynomials i?^ 1,K2 ' K3 (n(n+l)i) as well as that of the normalized spectral polynomial 
R n (n(n + l)t) converge to the measure p,Q defined in Theorem]]^ see Figure\J\ below. 




Figure 1 . Comparison of the theoretical density with the numer- 
ical density of the measure p,Q for Q(z) — z 3 — z. 

Remark 1. Note that equation ([3|) is usually lifted to the elliptic curve "J/ defined 
by y 2 = Q( z ) an d that on if? this equation takes the so-called Weierstrassian form 

- [n(n + l)p(u) + E]j S(u) = 0, 

where u is the canonical coordinate on the universal covering of ^ given by 

f 00 dw 

If Q(z) has all real roots then n(n + l)p(u) becomes a finite gap real periodic 
potential on WL, see [11]. Depending on the parity on n four of the eight spectral 
polynomials i?^ 1 '" 2 '" 3 (i?) determine the ends of bands of the spectrum of the cor- 
responding Schrodinger operator on M. Hopefully our results will find applications 
to the spectral theory of finite gap potentials. 

To end this introduction let us mention that in the present context one can actu- 
ally "guess" the last two formulas in Theorem [T] (iv) by WKB-type considerations, 
as we were kindly informed by K. Takemura in the final stages of this work. How- 
ever, such arguments fail to provide an accurate mathematical proof. By contrast, 
our methods use solely rigorous results involving orthogonal polynomials, elliptic 
integrals and hypergeometric functions, and we are currently unaware of any pos- 
sible "shortcuts" in this set-up. Finally, we should also note that so far neither our 
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methods nor WKB-type "guesses" could provide an answer to these questions in 
more general situations, see §3 below. 

Acknowledgements. We are grateful to M. Tater from the Institute of Nuclear 
Physics, Czech Academy of Sciences, for interesting discussions on this subject 
during his visit to Stockholm in November 2006. We would especially like to thank 
K. Takemura from Yokohama City University for numerous illuminating discussions 
prior and during his visit to Stockholm in September 2007. We would also like to 
thank the anonymous referee for his insightful comments and suggestions. 

2. Proofs 

Proof of Theorem^ Theorem [T] essentially follows from the main result of [7] after 
some amount of work. First we express the polynomial Sp n (t) as the characteristic 
polynomial of a certain matrix. In order to make this matrix tridiagonal so as to 
simplify the calculations we assume wlog that 

Q(z) = (z — e3)z(z — ei) = z 3 + vz 2 + wz 

with es < e% = < e\. Set 



T=(z 3 



d 2 d 

+ (az 2 + (3z + 7) — - 9 n (z - t), 



where v, w, a, (3, 7 are fixed constants and 9 ni t are variables. Assuming that S{z) 



a z' 



a\z 



n-l 



. . + a n with undetermined coefficients a,-, < i < n, we are 
looking for the values of 6 n ,t and aj, < i < n, such that T(S(z)) = 0. Note that 
T(S(z)) is in general a polynomial of degree n + 1 whose leading coefficient equals 
ao[n(n — 1) + an — 9 n ]. To get a non-trivial solution we therefore set 

n = n(n — 1 + a). 

Straightforward computations show that the coefficients of the successive powers 
z n , z n ~ x , . . . , z° in T(S(z)) can be expressed in the form of a matrix product M n A, 



where A — (ao, a\, 
matrix 



, a n ) and M n is the following tridiagonal (n + 1) x (n + 1) 



In, 2 





«n,2 
t-£n,2 
Tn,3 







7n,n+l 







with 



7rv 



v(n — i)(n — i + 1) + (3(n — i + 1) 
(n — i)(n — i + 1) + a(n — i + 1) 



ie {l,. 



1}, 
H}, 



(11) 



w(n - i + l)(n - i + 2) + 7(71 - i + 2) 



e{2,. 



!}• 



A similar matrix can be found in [10] and also in [19] . The matrix M n depends 
linearly on the indeterminate t which appears only on its main diagonal. If the 
linear system M n A — is to have a nontrivial solution A — (ao, a%, a n ) T the 
determinant of M n has to vanish. This gives the polynomial equation 

Sp n (t) = det(M n ) = 0. 
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The sequence of polynomials {Sp n (t)} n £i, + does not seem to satisfy any rea- 
sonable recurrence relation. In order to overcome this difficulty and to be able to 
use the technique of 3-term recurrence relations with variable coefficients (which 
is applicable since M n is tridiagonal) we extend the above polynomial sequence by 
introducing an additional parameter. Namely, define 

Sp n>i (t) = det M Hi i, i£ {1, ...,n + l}, 

where M n ^ is the upper i x i principal submatrix of M n . One can easily check (see, 
e.g., p. 20]) that the following 3-term relation holds 

SPn,i(t) = (t - £,n,i)Sp n ,t-i(t) - 1p n ,iSPn,i-2(t) , « G {1, . . . , n + 1}, (12) 

where £, ni i is as in (fTTj) and 

tpn,i = an.ijn.i, i€ {2, ...,«+ 1}. (13) 

Here we used the (standard) initial conditions Sp ny o(t) — 1, Sp n ^\(t) = 0. It is 
well-known that if all £ nj i's are real and all ip n /s are positive then the polynomials 
Sp n ,i(t), i £ {0, . . . , n + 1}, form a (finite) sequence of orthogonal polynomials. In 
particular, all their roots are real. Under the assumptions of Theorem Q] the reality 
of £ n ,i is obvious. Assuming that the positivity of ip n ,i is also settled (see Lemma Q] 
below) let us complete the proof of Theorem [T] For this we invoke [3 Theorem 1.4] 
which translated in our notation claims that if there exist two continuous functions 
£(t) and ij}{r), r G [0, 1], such that 

lim &,n=£(r), lim = #r), ^(r) > Vt G [0,1], 

i/(n+l)—fr >r 

then the density of the asymptotic root-counting measure of the polynomial se- 
quence {Sp n {t)} neZ+ = {Sp„, n +i(t)}„ eZ+ is given by 

l 



where for any x < y one has 

— , 1 if s £ \x,y], 

otherwise. 

From the explicit formulas for and ip n< i (see (fTT|) and (|13|) one easily gets 

C(r)= hm = - v (l - r) 2 , 

%l (n+l)—>T 

iP(r)= lim ^,„ = - W (1-(1-t) 2 )(1-t) 2 . 

i/(n+l)-»r 

Notice that the above limits are independent of the coefficients a, f3, 7 and that by 
the assumption on Q(z) made at the beginning of this section one has w — e^e\ < 0, 
which in its turn implies that ip( T ) > for r £ [0, 1]. The required density Pq{s) is 
therefore given by 



pq(s) 



dT 



dT 



TryJ-iw (1 - (1 - t) 2 ) (1 - r) 2 - (v(l - r) 2 + s) 5 
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where \f + is meant to remind that the integrand vanishes whenever the expression 
under the square root becomes negative. Introducing ^/v = 1 — r we get 



1 f 1 

pq{s) = ^L 



dv 



\J v\(Aw — v 2 )i> 2 — (Aw + 2vs)v — s 2 } 
In order to get rid of ^f~- + we rewrite 



Pq(s) 



1 



.(») 



dv 



2ttVv 2 - Aw Ju min (s) \Jv(v - v min {s))(v max (s) - v) ' 



(14) 



where v m in{s) and v m ax(s) are the minimal and maximal roots of of the equation 

(Aw - v 2 )v 2 - (Aw + 2vs)v - s 2 = 
with respect to v, that is, 

{v m in{s) = -(v 2 - 4w)~ 1 [2u; + vs + 2y/w{w + vs + s 2 )}, 

\v m ax(s) — -(v 2 - Aw)~ 1 [2w + VS — 2y / w(w + VS + S 2 )]. 

A few remarks are in order at this stage. First, by the real-rootedness of Q(z) 
one has v 2 — Aw — (e\ — e 3 ) 2 > 0. Second, we claim that 



£(t) - 2vW). £to + 2V^)J C [e 3 , ei] 

for t e [0, 1]. Indeed, 

ei - £(r) - 2V?W = ei + v(l - r) 2 - 2^/-w(l - (1 - r) 2 )(l - r) 2 



= ex(l - (1 - r) 2 ) - 2Vei(l-(l-r)2)V-e 3 (l-r)2 - e 3 (l - r) 2 



Vei(l-(l-r)2)- ^-^(l-r) 2 



> 



and similarly 



£(r) - 2v^(7) - e 3 = -«(1 - r) 2 - 2^(1 - (1 - r) 2 )(l - T ) 2 - e 3 

= ei (l - r) 2 - 2v/ ei (l-r) 2 v/- e3 (l-(l-r) 2 ) - e 3 (l - (1 - r) 2 ) 



Vei(l-r) 2 -v/- e3 (l-(l-r) 2 ) 



> 0. 



It follows that whenever s is such that the integrand in the above formulas is non- 
vanishing one has e 3 < s < e\ hence 

w(w + vs + s 2 ) = — e 3 ei(s — e 3 )(ei — s) > and 

2w + vs = ei(e 3 — s) + e 3 (ei — s) < 

for all s as above. Therefore 

v m in{s) + Vmax(s) = -2(2w + vs)(v 2 - Aw) > and 

Vmin(s)Vmax(s) = S 2 (v 2 - Aw) > 

from which we conclude that v max (s) > v m i n (s) > for e 3 < s < e±. 

Recall the definition and properties of the function uj from ((HI)-© and note that 



w(ei,e 3 ,s) 



(ei - e 3 )|s| 



|s|Vw 2 — Aw 



(15) 



ei(s — e 3 ) — e 3 (ei — s) 2w + vs 

so that v m in{s) and v m ax(s) may actually be rewritten as follows: 

|*| (l- v/l-^ei.es.s) 2 ) |s|(l + v/l-^(ei,e 3 ,s) 2 ) 



(ei - e 3 )w(ei,e 3 ,s) 



(ei - e 3 )w(ei,e 3 ,s) 
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Using these expressions combined with the fact that for any x < y one has 
n dv 2 




(which readily follows from 0) we deduce from (fH|) that 




e 3 )|s|u;(ei,e 3 ,s) \ ]j I - y/l - w(ei, e 3 , sf 
Now from ©-([7]) and the well-known identities T(l) = 1, Y{\) = \[k one easily 



gets 




v /l-^(e 1 ,e 3 , S ) 2 ^/l 1 ^ 2^1~cu(e 1 ,e 3 , S ) 2 \ ^ 



- e 3 )|s|w(ei,e 3 ,s) ^2 2 ' 1 - y/l - w(ei,e 3 ,s) 2 

by (fTT)|) . It is a remarkable fact due to Kummer and Goursat that for special choices 
of the numbers a, b, c the hypergeometric series F(a, b, c; z) obeys certain quadratic 
transformation laws. To complete the proof we need precisely such transformation 
properties, namely formulas 15.3.19-20 in [TJ with a = j, b = h and c = 1: 



The above identity for z = 1 — u>{e\, e 3 , s) 2 together with p7|) then yields 



/ \ / uj(e 1 ,e 3 ,s) /l 1 1 -u(ei,e 3 ,s) s . 

which by ([15")) amounts to 



y/2 [(ei + e 3 )s - 2eie 3 + (ii - e 3 )|s|] 

p/1 1 (ei +e 3 )s- 2eie 3 - (gi - e 3 )|s| \ 
1,2' 2' ' ( ei + e3 ) s -2e 1 e 3 + (e 1 -e 3 )| S |y' ■ 1 yj 

To prove the last two formulas of Theorem Q](iv) we use (|19[) and (J7J) with a = 6 = | 
and c = 1 in order to get an integral representation of Pq(s) in which one makes 
the following variable substitutions: 

e 3 (ei-x) ei(x-e 3 ) 
t = £ if e 3 < s < 0, t = if < s < ei. 

ei(x-e 3 ) e 3 (ei-x) 

The desired expressions are then obtained by straightforward computations. 

We thus established all formulas stated in Theorem [T] in the special case when 
e 3 < £2 = < e±. The general case reduces to this one by noticing that if e 3 < C2 < 
ei and Q(z) — (z — e 3 )(z — e-i){z — e\) then all the above arguments may be used 
for the polynomial Q(z + e^)- Hence the expressions for pq(s) given in Theorem[T] 
in the general case are obtained from (TT5)) - (flT)|) simply by replacing e\, e 3 and s 
with ei — 62, e 3 — C2 and s — respectively. □ 

Lemma 1. Let Q{z) — z 3 + vz 2 + wz and P(z) = az 2 + (3z + 7 be two polynomials 
such that Q{z) has three real distinct roots e\ > 62 = > e 3 and 

P(z) _ Cg C(2 a 3 

Q(z) z — e\ z — e 2 z — e 3 
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with ai, 02,0:3 > 0. Then the coefficients ^p n ^ defined in (|13p are all positive. 

Proof. One immediately gets a > 0, which implies that 9 n > for any positive 
integer n. Recall from (fTTj) and (fT3"|) that ip n ^ — a n ^ 7ly i. Now 

ct n ,i = 0^ 1 {i( n ~ — i + 1) — n(n — 1)] + a[(n — i + 1) — n]} 

and since both a and are positive it follows that a n ^ < 0. Let us show that under 
our assumptions one also has 7,^ < 0. For this it is clearly enough to show that 
both w and 7 are negative. Obviously, w — e^e\ < and since ct\, a 2 , 03 > the 
quadratic polynomial P(z) has two real roots interlacing with e3,e2,ei. Therefore 
P{z) has one positive and one negative root and positive leading coefficient hence 
7 = P(0) < 0. □ 

Proof of Theorem^ In order to deduce the differential equation satisfied by Pq(s) 
we note first that the restrictions of pq(s) to {e^,e2j and (e2,ei), respectively, are 
two branches of the same multi-valued analytic function (this can be seen e.g. from 
the last two expressions in Theorem [T]). Thus it suffices to derive the linear differ- 
ential equation satisfied by pq{s) restricted to, say, (e2,ei). Specializing formula 
(iv) of Theorem [T] to this case we get 



pq( s ) 



where 



2ny/(ei - e 2 )(s - e 3 ) 
dw 

"(i-«o(i- fc:air> 



and s € (e2,ei). By ([7]) we see that up to a constant factor Iq(s) is the same as 
the hypergeometric series F |, 1; |f^efyff^fj) ' which IS known to satisfy the 
following Riemann differential equation (see, e.g., [TJ §15.6]): 

T „, \ 1 ( 1 1 \ ri fo s , ( e 3 ~ e 2 )(e 3 - e{) 




Substituting Iq(s) = K^fs — e^pQ\s), where K := 2ny / ei — e 2 is but a constant, 
in the latter equation we get after some straightforward algebraic manipulations 

p' Q (s)+ ^~ Elie \ PQ ( S )^0, (20) 

which after multiplication by 8Q(s) coincides with the required equation (|10p . To 
calculate the indices recall from e.g. [T^] that for a second order linear differential 
equation p"{s) + ai(s)p'(s) + a,2(s)p(s) = its indicial equation at a finite regular 
or regular singular point § has the form 

C(C-l) + aiC + «2 = 0, 

where ct\ = lim s _,s(s — s)ai(s) and a 2 = lim s _>g(s — s) 2 a 2 (s). Thus for ([2H)l the 
indicial equation at each e*, i = 1, 2, 3, has the form C, 2 = 0, which implies that both 
corresponding indices vanish and that a solution to (|20|) might have a logarithmic 
singularity at any of these points, see Remark [2] below. The indicial equation of 
p"(s) + ai(s)p'(s) + a,2(s)p(s) — at 00 has the form 

C(C + 1) - fiiC + OL2 = 0, 

where a\ = lim^oo sai(s) and a 2 — lim s _,oo s 2 a 2 (s). Therefore the indicial equa- 
tion at 00 for (f2T)]) is ((( + 1) — 3£ + § = whose roots are \ and §. □ 
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Remark 2. It is not difficult to show that 

1„„/ 16(ei-e 2 )(e 2 -e 3 ) ' 
1 l0 &{ ( ei _ e3 )| s -e 2 | 

Pq( s ) ~ 7T- /7 f7 . as s -» e a , 

2?r V( e i - e 2)(e2 - e 3 ) 

so that pq{s) always has a logarithmic singularity at e%. 

Proof of Corollary [H The main idea of the proof of Corollary [T] is that the poly- 
nomial part of any Lame solution to ((3]) itself satisfies a very similar differential 
equation. Let us illustrate this in the case of Lame solutions of type 2. (The other 
cases can be dealt with in the same way.) Assume that S(z) is a polynomial of 
degree n such that S(z) := (z — ei) 2 (z — e 2 ) T 5(z) solves the equation 



4Q(z) 




S'(z) 



+ V(z)S(z) = 



for some linear polynomial V(z). Recall that the (unique) root of each such V(z) 

- - o 

is also a root of the scaled spectral polynomial Rn 2 ' (n{n + l)t). Substituting 
S(z) — (z ~ ei) z (z — 62)2 S(z) one gets after some straightforward calculations that 
S(z) satisfies the equation 



AQ(z) 



S"(z) + \ f_L_ + _L_ + _L_ )5'(z) 

2 \2— e\ z — e2 2: — e3 



+ F(z)S(z) = 0, (21) 



where V"(z) = V(z) + (z — ei) — (z — 62) + 4(z — 63). Notice that in the coefficient 
in front of S'(z) the numerators of the first two simple fractions are increased by 
1 while the remaining one is unchanged. Thus we are practically in the situation 
covered by Theorem [TJ The only (slight) difference is that we are not considering 
the asymptotic distribution of the roots of polynomials V(z) but that of V(z). 
However, since V(z) — V(z) + (z — e{) — (z — e 2 ) + 4(z — e 3 ) and the leading 
coefficient 6 n of V(z) tends to 00 as n — > 00 it follows that both these families 
of polynomials actually have the same asymptotic root distribution. Finally, note 
that by (jlj the scaled spectral polynomial R n (n(n + i)t) equals the product of 
four corresponding polynomials R* 1,K2 ' K3 (n(n + l)t) and that all eight polynomial 
families {i?^ 1 '" 2 ' 143 (n(n + l)i)} ng z + have the same asymptotic root distribution 
when n — > 00. This implies that the family {R n (n(n + l)t)}„ e z + itself has the 
same limiting root distribution. □ 

3. Remarks and Conjectures 

1. So far we were unable to extend our method of proving Theorem[T]to the case of 
a cubic polynomial Q(z) with complex roots. The main difficulty comes from the 
fact that Theorem 1.4 of [7] seems to fail in this case. Indeed, if it were true then the 
support of the resulting root-counting measure would be two-dimensional. However, 
numerical experiments strongly suggest that this support is one-dimensional, see 
Figured 

Problem 1. Generalize Theorem[T]to the case of a complex cubic polynomial Q(z). 

Based on our previous experience (comp. [4J, [5] ) we conjecture that: 

(i) The asymptotic root distribution of the Van Vleck polynomials for equation 
(f5|) with complex coefficients is independent of the a^s and depends only 
on the leading polynomial Q(z); 

(ii) Denoting the above limiting distribution by [iq we claim that its support 

is straightened out in the canonical local coordinate w(z) = V~ , . 

Jz ° \/Q{t) 



12 



J. BORCEA AND B. SHAPIRO 



1[* - 
0.8 
0.6 
0.4 
0.2 

Ot , , L. , , , ' ' ' ■■ : 

-0.4 -0.2 0.2 0.4 0.6 0.8 1 

Figure 2 . The roots of the scaled spectral polynomial Sp^o (t) for 
Q(z) = z(z-l)(z + i-i). 

2. Even more difficulties occur when dealing the more general Lame equation ^ 
since in this case Van Vleck polynomials and their roots can not be found by 
means of a determinantal equation. They arc in fact related to a more complicated 
situation when the rank of a certain non-square matrix is less than the maximal 
one, see [5]. 

Problem 2. Describe the asymptotic distribution of the roots of all Van Vleck 
polynomials for equation ([2]) when n — > oo. 

An illustration of this asymptotic distribution is given in the next picture. 




Figure 3. The union of the roots of 861 quadratic Van Vleck 
polynomials corresponding to Stieltjes polynomials of degree 40 for 
the classical Lame equation Q(z)S"(z) + ^^-S l {z) + V(z)S(z) = 
with Q(z) = (z 2 + l)(z - 3i - 2)(z + 2i - 3). 
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